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We model the quantum Hall effect in heterostructures made of two gapped graphene stripes with 
different gaps, Ai and A2. We consider two main situations, Ai = 0, A2 7^ and Ai = — A2. They 
are different in a fundamental aspect: only the latter feature kink states that, when intervalley 
coupling is absent, are protected against backscattering. We compute the two terminal conductance 
of heterostructures with channel length up to 430 nm, in two transport configurations, parallel 
and perpendicular to the interface. By studying the effect of disorder on the transport along the 
boundary, we quantify the robustness of kink states with respect to backscattering. Transport 
perpendicular to the boundary shows how interface states open a backscattering channel for the 
conducting edge states, spoiling the perfect conductance quantization featured by the homogeneously 
gapped graphene Hall bars. Our results can be relevant for the study of graphene deposited on 
hexagonal Boron-Nitride as well as to model graphene with an interaction-driven gapped phase 
with two equivalent phases separated by a domain wall. 

PACS numbers: 73.43.-f, 73.22.Pr, 72.80.Vp 
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I. INTRODUCTION 

The Hall conductance in quantum Hall bar is so ac- 

2 
curately described by (Jxy — rij^, where n is an integer 

number, that it is used^ as our standard definition of 
the ratio of such fundamental constants as the square of 
the electron charge e^ and the Planck constant h. The 
origin of this extraordinary quantization, by which the 
conduction properties of a device are independent of the 
material properties, is intimately linked to the fact that 
in quantum Hall bars transport takes place only through 
the edges which host chiral states for which backscatter- 
ing is forbidden.EEl In turn, the existence of chiral edge 
states that permit non-dissipative transport is warranted 
by the topological order of the electronic states of the 
two dimensional gas states. For two decades, this state 
of affairs was only observed at sub-Kelvin temperatures, 
under high applied magnetic fields, in two dimensional 
electron gases, hosted by carefully designed modulation 
doped semiconductor heterostructures. The discovery of 
quantum Hall effect on graphene^ on one side, and the 
proposal^HSl and subsequent discoverjE! of quantum spin 
Hall insulators on the other, have dramatically expanded 
the materials and experimental conditions under which 
non-dissipative quantum transport linked to topological 
order can occur. 

Most of the striking electronic properties of graphene 
are related to the absence of a gap separating the conduc- 
tion and valence bands which can t hereby be described 
in terms of massless Dirac fermions.l^SIll] in particular, 
the magneto-electronic properties of graphene are funda- 
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mentally different from a non-relativistic two dimensional 
electron gas on three counts; ^^ the existence of two iden- 
tical sets of Landau levels, for electrons and holes, the 
scaling of their energy with VB, as opposed to linear 
scaling of non-relativistic fermions, and the existence of 
the n — Landau level with zero energy. T hese proper- 
ties make the quantum Hall effect in graphene^El different 
from the one originally discovered in GaAs two dimen- 
sional electron gases.'^ 

There are several physical scenarios that motivate the 
study of the electronic properties of gapped graphene. 
First, a gap could be opened by interaction driven elec- 
tronic order ,11211151 specially when a high magnetic field is 
applied. Seco nd, as a result of the influence of the sub- 
strate like SidiSliZlor hexagonal Boron Nitride (BN),'l21 
although the lattice mismatch is known to complicate 
this second possibility.lSHlIl Third, BN itself can be de- 
scribed with the tight-binding Hamiltonian of gapped 
graphene and the k ■ p Hamiltonian of other two dimen- 
sional materials with hexagonal symmetry, such as M0S2, 
can be described with a massive Dirac Hamiltonian.'^SIlll 
Fourth, i ntri nsic spin orbit coupling also opens a gap in 
graphene,!^ with different sign at the two valleys, albeit 
very smalLl^Sl 

All of this leads to the question of how magnetotrans- 
port properties of graphene and graphene-like materials 
change when a gap opens or, in the long wavelength limit, 
how massive and massless Dirac fermions are different in 
their reaction to an applied magnetic field. It turns out 
that, when the gap is opened by a constant staggered po- 
tential, i.e., a potential that acts with opposite sign in the 
two sublattices of graphene, as it happens if pseudospin 
magnetism, or in the case of hexagonal Boron Nitride, 
the answer to the question is quite straightforward from 
the theory standpoint. This occurs because there is a 
simple one-to-one relation between the energy levels and 
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FIG. 1; (Color online) Scheme of the two geometries consid- 
ered in the text: (a) perpendicular transport to the interface 
between the two different gap regions and (b) for the parallel 
transport. Arrows indicate edge, interface and kink states in 
the quantum Hall regime, (c) Detail of the armchair ribbon 
Nw ^ 7 {W ^ 0.74 nm) and iVi = 3 (L = 1.3 nm) and 
(d) zigzag ribbon edges Nw = 4 {W — 0.9 nm) and Nl ~ 5 
(L = 1.2 nm). 



wave functions of a bipartite lattice Hamiltonian with no 
staggered potential and those of the s ame l attice when a 
constant staggered potential is added.'^SlHZl This relation 
permits to anticipate that quantum Hall effect of massive 
Dirac fermions is much closer to the one of massless Dirac 
fermions than to the one of non-relativistic electrons.'23 

The situation becomes more interesting when the gap 
-or mass-, is not homogeneous. This could be the case, 
for instance, of a heterostructure made of two graphene 
or graphene-like materials with different gaps Ai and A2, 
such as the atomic layers of hybridized BN and graphene 
domains,!^ or if the gap is substrate-driven and, due to 
lattice mismatch, feat ures a mplitude modulations larger 
the graphene unit cell.^^ ^^ Here we study the electronic 
properties of heterojunctions formed between two gapped 
graphene-like systems, with different gaps, in the quan- 
tum Hall regime. When decoupled, both Hall bars would 
have its own set of chiral edge states. When coupled, the 
way edge states merge determines the electronic proper- 
ties of the interface, which is the focus of this work. 

As we discuss below, we encounter two types of inter- 
face states. The hybridization of pre-existing edge states 
that carry electrons in opposite directions will give rise 



to interface states that can carry electrons in both direc- 
tions. In contrast, the merger of two bars with opposite 
gaps gives rise to two interface states whose energy lies 
within the gap. At zero magnetic field these states can 
be rationalized as follows. When restricted to one valley 
T, graphene electrons can be assigned a Chern number 



2A 



, with T = ±1. According to the index theorenP^l^ 
the interface between two insulators with Chern num- 
bers Til and n2 should host Q — \ni — n2\ domain wall 
states, which will correspond to zero modes in the case 
of one dimension--'^ and chiral bands or kink states in 
two dimensions. These are very similar to the recently 
discussed kink states in the interface of two graphene 
bilayers^^HSU ^nd multilayers^ with a gap opened by the 
application of an electric field with opposite direction at 
the two sides of the junction. 

The rest of this paper is organized as follows. In 
section II we review the electronic structure of gapped 
graphene under the influence of a perpendicular mag- 
netic field, both within the conventional tight-binding 
model as well as its long wavelength counterpart, the 
massive Dirac fermions. In section HI we study the elec- 
tronic structure of graphene heterojunctions in a geome- 
try that preserves translational invariance along one di- 
rection, which simplifies the discussion and permits to 
unveil the appearance of interface states. In section IV 
we study the quantum transport properties of the kink 
states in these structures, including the effect of disorder. 
In section V we consider heterojunctions formed by two 
semi-infinite Hall bars made of gapped graphene-like sys- 
tems with different gaps. Transport in this type of het- 
erojunction could be used to probe the interface states, 
which enables backscattering between the otherwise de- 
coupled chiral edges states. In section VI we summarize 
our conclusions. 



II. STRIPES OF GAPPED GRAPHENE IN A 
MAGNETIC FIELD 

A. Tight Binding Hamiltonian 

Here we review the well studied tight-binding model for 
grap hene unde r the influence of a perpendicular magnetic 
ggjjjnnilSZllll ^g consider the standard tight-binding 
model for graphene, with one orbital per atom and first 
neighbor hopping t sa 2.7 eV, with a staggered poten- 
tial A{R) that might be position dependent. A given 
graphene unit cell, located at R, has two atoms, denoted 
by A and B. Using this notation, the zero field tight- 
binding Hamiltonian reads: 



H 
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where t^ ^, is nonzero only for first neighbors and a^, 6^ 
annihilate an electron at the A and B sites of unit cell R 



defined on a honeycomb lattice, and taking A as a con- 
stant along the entire system, this Hamiltonian describes 
graphene with a gap of 2 A in both valleys. In the rest of 
this paper the spin degree of freedom is ignored. Results 
for non-interacting electrons with spin can be obtained 
by adding the Zeeman shift to the obtained bands. 

Within this tight-binding description, the effect of the 
applied magnetic field is included replacing the hopping 
ti.2 between sites 1 and 2 of the lattice of the B = 
Hamiltonian by ii^2 -^ ii,2e'*i-^ wherePHal 



*1,2 = 



A- dr, 



(2) 



is the circulation of the vector potential A associated to 
the magnetic field B and the labels 1 and 2 stand for 
the coordinates of the two atoms whose hopping integral 
is being calculated. This is the lattice analogous of the 
canonical substitution for the free electrons, where the 
momentum operator pis replaced by p— eA. Notice that 
the phase $1.2 that modulates the hopping is propor- 
tional to the ratio of the magnetic flux per unit cell and 
the magnetic flux quantum $0 = -. 

In the following we assume that graphene lies in the 
z — plane, and we take B = B{0, 0, 1). Taking advan- 
tage of the gauge symmetry, we choose 



so that 



$1.: 



A = B(-y,0,0), 
eB (x2 ~xi)(j/2 +Vi) 

h 2 



(3) 



(4) 



where {xi,yi) are the cartesian coordinates to atoms 1 
and 2. With this choice, the Hamiltonian keeps transla- 
tional invariance along the x direction. 



B. Effective mass approximation 

Whereas the tight-binding approach provides a fairly 
complete description of the non-interacting electrons in 
graphene under the effect of a magnetic field, as we dis- 
cuss below, most of the results for states with energies in 
the neighborhood of the Dirac points can be rationalized 
by making use of the k ■ p description of the bands.!^ 
The magnetic field introduces a new length scale in the 
problem: 



h 



(5) 



We assume a sufficiently high magnetic field so that Ib < 
W, where W is the width of the ribbon, and hence the 
bulk quantum states become localized, the spectrum of 
states away from the edges becomes discrete, the bulk 
is an insulator, and dispersive and conducting states are 
only possible at edges. For typical magnetic fields, we 



also have a « Ib, where a = ^/3acc is the graphene 
lattice constant, which enable a description of the energy 
levels in terms of an effective k ■ p Hamiltonian. 

The effective k ■ p or effective mass Hamiltonian turns 
out to be isomorphic to the Dirac Hamiltonian at the two 
valleysP 



Hr = vp {HxCx + rllya 



y'-'v) 



Ac 



(6) 



where li = p — eA is the canonical momentum opera- 
tor, Vp = 'itacc/'^^, CT are the Pauli matrices describing 
the graphene sublattice degree of freedom and t = ±1 
describes the valley index. 

Using the gauge defined in Eq. ^ leads to 



Hr 



A Vp [Px - eBy 

Vp [px - eBy - irpy] -A 



iTpyl 



,(7) 



This Hamiltonian is translationally invariant along the x 
direction, so that we can assume its eigenfunctions are 
products e^''''^^(j)n{kx,y) which permit replacing the op- 
erator Px by the quantum number hkx in Eq. ([7]). We 
define the dimensionless canonical operators: 



Qikx) 



kxi'b I 5 



and 






(8) 



(9) 



combined with the intrinsic energy scale associated to the 
Fermi velocity. 



TiujQ _ hvp 



(10) 



Notice that the role of kx is to shift the eigenvalues of 
the Q{kx) operator. It is very convenient to define the 
ladder operators: 



a{kx) = \[Q{kx)+iP), 
v2 



(11) 



which satisfy [a{kx) , a{kxy^] = 1. For simplicity, in the 
following we omit the kx dependence of the a operator. 
We thus can write the Hamiltonian for the r = -t-1 valley 
as: 



H+ikx,Py) 



A ^%ackt 

\/2 



V2 



-A 



(12) 



whereas for the r = — 1 valley the Hamiltonian reads 

H^{kx,Py) = H{kx,Py)'^. 

In order to find the eigenfunctions and eigenvalues of 
these Hamiltonians, it is convenient to compute their 
square: 



Hr=±l - 



H^ir) 
n^-T) 



(13) 



where H^{t) = A^ + i {Tu^of {a'^ a + ^). 

In the following we denote the the eigenstates of the op- 
erator a^a as 0„, with eigenvalues n = 0, 1, ... and eigen- 
functions (/)„. The states 4'n are the standard harmonic 
oscillator wave functions, centered around y — kxls- The 
eigenstates of H'^^^.i, and thereby eigenstates of "Hr, de- 
noted by (pn, fall in two categories: the so called, zero 
Landau level, with a sub-lattice polarized wave function, 
and the normal Landau levels. 



Landau Levels 



From Eq. ( 13 ) it can be seen right away that eigen- 



states can be written as: 



Ki^) 



Ar. 



B 



n'rii — T 



(14) 



where A„ and i3„ are coefficients that are determined 
by requesting that (j)n are also eigenstates of the Dirac 
equation. The corresponding eigenenergies are Ef^ — 
A^-t-i {huJo) {n + ^^)- Therefore, the general equation 
for the eigenvalues of the Dirac Hamiltonian under t he in- 
fluence of a perpendicular two dimensional field areP^^Sl 



£;„(r) = ±JA2 



2 (^^o)^ 



1-r 
2 



(15) 



It is apparent that the Landau level energies are inde- 
pendent of kx- Therefore, they give rise to flat bands with 
a very large degeneracy. Moreover, there is an additional 
twofold valley degeneracy given by: 



^„(r = -l) = S„+i(r = +l) 



(16) 



Of course, Eq. (14) is only mathematically defined if 
both n and n + t are non negative. As a results, the 
minimal value that n can take is n = 1 for t = —1, 
and n = for r = -f 1. For these states, the energy can 

be written as -Ejv — ^V^^ + 5 i^^o) [^ + l)i where 
A^ -(- 1 is a strictly positive integer. In summary, these 
states come in doublets, on account of the valley degree 
of freedom, and in addition have electron-hole symmetry. 



2. Zero Landau Level 

In addition to these states, for each valley there is 
an extra eigenstate of H^ with the minimal eigenvalue 
H^Zt- = A^iV- They are: 



and 



Zt=+1 — 



Zt=-1 = 
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FIG. 2: (Color online) (a) Scheme of the Landau levels for 
gapped graphene in the fc-p approximation. Notice the valley 
dependence for the n — Landau level. Labels show the wave 
function resolved in sub lattice component, (b) Scheme of the 
Landau levels in graphene ribbon, where valley mixing takes 
place. The energies of the levels correspond to A = 0.02 t 
and $ = 0.0002 [B = 15.8 T). 



We thus see that these wave functions are very special: 
they are sublattice polarized. It can be verified right 
away that these wave functions satisfy: 



Ht-Zt- = tAzt-. 



(19) 



Thus, the energy of the zeroth Landau levels becomes 
valley dependent due to the mass term A, as shown in 
Fig. 2(a). We can relate this to the fact that the mass 
introduces an orbital magnetic moment with valley de- 
pendent orientation.HD 

A summary of the energy spectrum for gapped 
graphene electrons under the influence of a perpendicular 
magnetic fleld, described within the fc-p approximation is 
shown in Fig. [2J Whereas for all the flnite \n\ levels each 
valley contributes with one Landau level, so that they 
come in couples in graphene, the n — Landau levels are 
valley polarized, so that there is only one for the electron 
sector and one for the hole sector. For A = these two 
n = levels are degenerate. However, this degeneracy is 
lifted for gapped graphene, and a gap between them is 
open. This discussion has ignored the spin degree of free- 
dom, which would add an additional twofold degeneracy 
to all the levels, broken by the Zeeman splitting. 



C. Landau levels and edge states 

We now consider the spectrum of the edge states of 
gapped graphene in the quantum Hall regime. Edge 
states are important in this regime because they provide 




states localized away from the edges. We have verified 
that the energies are described by Eq. ( 15 1 . In particular, 
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FIG. 3: (Color online) Band structure in the quantum Hall 
regime for uniform massless armchair (a) and zigzag (c) 
nanoribbons. And of the massive (A = 0.02 t) ribbon arm- 
chair (b) and zigzag (c). In the four cases the magnetic field 
is fixed $ = 0.0002 {B = 15.8 T) and the dimension are 
Nw = 1000 [W = 123 nm for armchah W = 213 nm for the 
zigzag). 



the only transport channel. Whereas it is possible to pro- 
vide an approximate description for edge states within 
the effective mass k ■ p approximation used for the bulk 
states in the previous section, here we apply the tight- 
binding methodology to compute the the energy levels 
of graphene stripes o f wid th W under the influence of a 
strong magnetic field.'^S^ 

The use of the gauge choice of Eq. ([3|, permits studing 
quantum Hall bars that are infinite along the x direction 
and have a finite width along the y axis. For a unit 
cell with N atoms, we obtain A'' bands en{kx)- We can 
consider two geometries, with either zigzag or armchair 
edges (see Fig. [ijc-d)).'*^ In the rest of the paper we de- 
scribe the magnetic field in terms of the magnetic flux per 

hexagon in the honeycomb lattice, <f> = 2$°^^ normal- 
ized to the magnetic flux quantum $0 = ■ For reference, 
a normalized magnetic flux of <& = 10~^ corresponds to 
B = 7.9 T. 

In Fig. p^we show the energy bands, denoted by e„(fc), 
for two different graphene stripes, with zigzag and arm- 
chair terminations, for $ = 2 x 10"^ and A = (Fig. 
[3];a) and (c)) and A = 0.02 t (Fig. [sfb) and (d)). The 
width of the ribbons is W = 213 nm for the zigzag and 
W = 123 nm for the armchair. There are several things 
to notice. Bands are flat in wide regions of the Brillouin 
zone and dispersive otherwise. An analysis of the corre- 
sponding wave- functions indicate that, except in one case 
described below, flat bands correspond to Landau level 



the energy gap between different Landau levels is index 
dependent, as expected for Dirac electrons and different 
from non-relativistic electrons. 

The dispersive states correspond to states localized at 
the edges. There is a linear relation between localiza- 
tion along the transverse direction of the ribbon and the 
momentum k^, as expected from Eq. (IS]). The edge 
velocity w„(fc) — ^ '^g^; ' changes from one edge to the 
other. The emergence of these chiral edge states whose 
energy lies in the gap between Landau levels anticipates 
the very peculiar quantized transport properties of the 
system, characteristic of the quantum Hall effect.^ 

In agreement with the effective mass results, the flat 
bands (Landau levels) have a twofold degeneracy, both 
in the armchair and zigzag cases, except for the n — 
level. In the case of armchair termination the degeneracy 
occurs at the same k point, whereas in the case of zigzag, 
there are two sets of bands, that can be ascribed to the 
two valleys.HS The A term shifts the position of all the 
Landau levels and splits the n — levels opening a trans- 
port gap even at the edges, also in agreement with the 
effective mass results. The presence of a two flat bands, 
at a given valley, associated to the n = Landau is in 
apparent conflict with the effective mass approximation 
(see Fig. [2]). It turns out that one of the two n = 
flat bands at each valley, is an edge state, rather than 
a bulk state.'^ This statement is further clarified in the 
next section. 



III. ELECTRONIC PROPERTIES OF 

GRAPHENE QUANTUM HALL BARS WITH 

INHOMOGENEOUS A 

We are now in position to study the electronic states 
of graphene quantum Hall bars where A is not constant. 
For that matter, we consider the simplest situation, a 
ribbon of width W where top and bottom halves have a 
different mass At and A^ . We consider two cases At = 
—As and A^ 7^ 0, A^ = 0, both armchair and zigzag 
terminations (four cases in total). The band structures, 
together with the local current density: 



Jxiy,E) 



E 



I0n,fejy)|' 



, oen(kx) 
dkx 



5{E~e^{kx)), (20) 



are shown in Fig. [4] 

For a given Landau level n with wave function ipnik, y), 
and within a given valley, there is a relation between 
the k quantum number and the average vertical position 
{'ipnik,y)\y\'ipnik,y))- Thus, plots of the current density 
as a function of y provide complementary information to 
bands tn{k). In panels (a-b) we show the armchair rib- 
bon with At 7^ 0, A^ = 0. It is apparent that, according 
to their location in the ribbon, we can distinguish three 
types of states: edge states, bulk states and, in contrast 
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FIG. 4: (Color online) Band structure (left row) and the cor- 
responding current density defined by Eq. (201 (right row). 



Here we consider transport parallel to the interface (see Fig. 
[l](b)) for armchair (a-d) and (e-h) zigzag ribbons. We fix the 
mass of the upper half of the ribbon At = 0.05 t, and the 
lower half of the ribbon has either Ab = ((a-b) and (e-f)) 
or Ab ~ —At (panels (c-d) and (g-h)). The magnetic flux 
is $ = 1.1 X 10"^ {B = 87 T) and the width of the ribbon 
is Nw = 300 {W = 36.9 nm for armchair W = 63.9 nm for 
zigzag) 



This hybridization results in two intertwined oscillating 
bands. For n = Landau levels, of the interface states 
are unique and join the two pairs of n = Landau lev- 
els. They can also be interpreted as regular edge states 
of the massless half confined by the gap on one side and 
the vacuum on the other. 

Results become more interesting for the armchair rib- 
bon with At = —Ab ^ 0, shown in Fig. Wc-d). Bulk 
Landau levels, edge states and n ^ interface states 
are very similar to the previous case. The main dif- 
ference occurs for the interface states for the n = 0, 
which fill the gap almost completely. The two counter- 
propagating interface states undergo a small anti-crossing 
at zero energy. These interface states that reside in the 
gap are quite similar to the kink states reported for bi- 
layer graphene with a position dependent off-plane elec- 
trie field.l^MIl 

The discussion for zigzag ribbons goes along the same 
line. For a zigzag ribbon with A^ ^ 0, Ab — 0. (Fig. 
W[c) and Fig. Illf)) we have two replicas of the Landau 
levels and their edge/interface states for each valley. For 
the n 7^ Landau levels there are bulk flat bands, and 
dispersive edge and interface states, very much like in the 
case of armchair ribbon. For the n — Landau levels. 
Fig. |4] shows four type of bulk states, attending to the 
sub-lattice (cr^) and valley (r^) indexes: i) gapless with 
{a^ = +l,T^ = +1), n) gapless with {a^ = -1,t^ = -1), 
both with zero energy. Hi) and iv) gapped (A > 0), with 
energy r^A and either (ct^ = +1,t^ = +1) or (cr^ — 
— l,Tz = —1), as expected from the effective mass theory. 
In addition, at each valley there is a kink state that joins 
the gapless n — Landau level, with the corresponding 
gapped state. This kink states share spectral range with 
edge states. 

The properties of the zigzag ribbon with At = — A^ 
are in line with the other cases. The main feature here is 
the presence of a kink state at each valley that, in con- 
trast with the armchair ribbon, it has no anti-crossings 
and covers the gap completely. The wave functions of 
the kink states are located at the interface, as expected. 
The velocities of the kink states are opposite for valley. 
This is one of the main results of this manuscript: we 
predict the existence of counter-propagating valley po- 
larized states at the interface of two graphene quantum 
Hall bars with opposite masses. 



with ribbons with homogeneous A, interface states lo- 
cated at the boundary between the massive and massless 
sectors. For the bulk states we obtain two different set 
of flat Landau levels, corresponding to the massless and 
massive halves, respectively. The edge states are quite 
similar to those of the homogeneous mass case. For the 
n y^ Landau levels, the interface states can be inter- 
preted as the hybridization of the two pairs of counter- 
propagating edge states from the homogeneous A sectors. 



IV. TRANSPORT PROPERTIES OF KINK 

STATES 

We now discuss the robustness with respect to disor- 
der of transport properties of the kink states found in 
ribbons with Ay = -Ab- For that matter, we consider 
the geometry shown in Fig. [l} an infinite stripe with 
three regions: two semi-infinite electrodes without disor- 
der joined by a central region, of length L, that features 



an Anderson disorder potential™ 



(21) 



where Vi is a random variable uniformly distributed over 
the interval [— Vb, Vb], where the energy scale Vb sets the 
strength of the disorder potential. The transmission is 
calculated for each disordered configuration and averaged 
over a different disorder configuration realizations. 

Making use of the partition method and the Green 
function approach, outlined in the appendix, we com- 
pute the scattering transmission function T[E), which re- 
lates to the two terminal elastic conductance through the 

2 

Landauer formula G — %-T{Ep). The transmission func- 
tion is the sum over of the transmission coefficients Tn of 
the channels n available at a given energy. For an ideal 
transmission channel without backscattering, T„ = 1. A 
completely blocked channel gives T„ = 0. 

The edge states in quantum Hall bars are the canonical 
example of ideally transmitting channels, with T„ = 1. 
This leads to a quantized two terminal conductance, 



G 



where n is an integer number. In the four 



terminal geometry the two terminal longitudinal con- 
ductance vanishes, but the Hall conductance quantiza- 
tion is related to the quantization of the two terminal 
conductance. ■ 

We consider first the transport properties of the arm- 
chair ribbon with opposite mass Ay — —A^. Since the 
calculation of the transmission coefficients requires the 
determination of the Green function of a system with 
2 X Nw X Nl atoms, and we consider lengths Nl up 
to several thousands, it is computationally convenient to 
choose a smaller Nw, but large enough so that there is no 
interedge coupling. This also makes it necessary to take 
unrealistically large values of B. However, we expect that 
the simulated structures have the same properties than 
wider ribbons with smaller magnetic fields. 

We study now the transmission as a function of energy, 
for a fixed length of the disordered region N^, shown 
in Fig. p[c). The shadow background with quantized 
steps shows the transmission without disorder, and the 
black line is the transmission for a length Nl = 30 for 
Vb = 0.1 t. The stepwise function for Vb = 0, i.e., in the 
absence of disorder, reflects the number of modes at a 
given energy, starting from when the energy lies inside 
the gap (coming from anticrossing of the kink states) and 
increasing as the Fermi energy reaches new edge states. 

The effect of disorder in Fig. ISlc), for a fixed channel 
length, is clearly energy dependent. This is more clearly 
seen in Fig. p[d), where we plot the conductance, aver- 
aged over 10 disorder configurations, at the four energies 
marked in Fig. [sja, c), as a function of the length of 
the ribbon. The vertical error bars reflect the standard 
deviation. It is apparent that, as the length of the trans- 
port channel L increases, the transmission tends to the 
quantized values, when the energy crosses edge states 
(green line), but tends to zero in the case of the kink 
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FIG. 5: (Color online) Electronic properties of an armchair 
ribbon with a stepwise constant gap At = —As = 0.3 t (see 
Fig. 0(b)), under a magnetic field $ = 0.0081 {B = 640 T), 
Nw = 100 {W = 12.3 nm). (a) Energy bands, (b) Current 
density map for the disorder-free structure, (c) Two terminal 
conductance as a function of Ef (in units of the hopping t) 
with disorder (Black line). As reference the conductance of 
the disorder-free ribbon has been included as a shadow region, 
(d) Two terminal conductance as a function of channel length 
for four different energies (marked in (a) and (c)). In the 
disorder cases (c-d) the Anderson parameter is Vb = 0.1 t, an 
average over 10 disorder configurations was performed and 
the error bars reflect the standard deviation. 



states (black and red lines). Furthermore, when the en- 
ergy crosses both edge and interphase states (blue line) 
the channels corresponding to the interphase states are 
attenuate as the length is increased. Therefore, backscat- 
tering is possible, which is expected since two states with 
opposite velocities coexist in the interface region. 

When the same analysis is done for a zigzag ribbon 
with Ax = — A^, the results found are quite similar 
(see Fig. [6]). Therefore, our transport results show that, 
unlike edge states, kink states living in the interface of 
two gapped graphene regions with opposite gaps are not 
protected. It must be stressed that, very much like in 
the case of gapped bilayer graphene, the existence of two 
valleys prevents the robustness of kink states with respect 
to disorder, as it would happen in the case of kink states 
located at the domain wall for a two dimensional electron 
gas of Dirac electrons. 
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FIG. 6: (Color online) (a), (b), (c) and (d) as in Fig. [s] for the 
case of a zigzag ribbon with a stepwise constant gap A with 
At = —As = 0.3 t, under a magnetic field with $ = 0.0081 
[B = 640 T), Nw = 100 {W = 21.3 nm). 



V. QUANTUM HALL TRANSPORT 

PROPERTIES OF GAPPED GRAPHENE 

HETEROJUNCTIONS 

In this section we study transport across of the bound- 
ary that separates two gapped graphene Quantum Hall 
bars. Our main goal is to analyze the backscattering 
that takes place at the interface due to interface states 
that connect the edges. A scheme of this heterostructure, 
made of two regions with constant mass A^ = — A/j = A 
that meet at an abrupt interface, is shown in Fig. RTa). 
We consider transport across the interface region which, 
in addition to the stepwise constant mass has Anderson 
disorder with Vo = 0.1 t over a distance Nl = 30. Away 
from the interface each side of the junction hosts a set of 
Landau levels and edge states as described in the section 
II. Because both electrodes have a gap 2A, conductance 
is only possible for states with energy _E > A. 

Both sides of the structure host chiral edge states that 
do not backscatter, even in the presence of disorder, un- 
less electrons at one edge can undergo scattering to the 
other edge. This could be enabled if interface states, i.e., 
states running along the interface, perpendicular to the 
transport direction, are available at the Fermi energy and 
are not completely blocked by disorder. Below we show 
that interface states, perpendicular to the transport di- 
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FIG. 7: (Color online) Electronic properties of a graphene 
heterostructure made of two semi-infinite armchair ribbons 
with opposite mass A^ — — A^ = 0.3 t (see Fig. IT] (a)), un- 
der a magnetic field $ = 0.0081 {B = 640 T). (a) Current 
density map corresponding to the zigzag infinite ribbon along 
the interface, (b) Two terminal conductance as a function of 
Ef with disorder (Black line) for a central region of length 
Nl = 30 (L = 12.8 nm) and width Nw = 50 (W = 6.1 nm) 
and average over 10 disorder configurations was performed. 
As reference we include the conductance of the disorder-free 
ribbon has been included as a Grey line for A_l = —An and 
the homogeneous mass Al = A_r as a Green shadow, (c) 
Two terminal conductance as a function of the ribbon width 
for four different energies (marked in (b)) calculated with dis- 
order, a channel length Nl = 30, averaged over 100 config- 
urations. All the disordered cases have been calculated with 
Anderson disorder Vb = 0.1 t. 



rection, can act as an efficient shortcut between the right- 
goers in one edge and the left-goers in the other (see Fig. 
[T]), providing a backscattering channel which destroys the 
conductance quantization. It must be stressed that, in 
this geometry, the in-gap kink states studied in the pre- 
vious section, do not play a role in transport, due to the 
absence of transport states in the electrodes inside the 
gap. 

In Figs. [7] and [8] we show results for the conductance 
of both armchair and zigzag heterojunctions. The in- 
terface between the two armchair (zigzag) ribbons with 
opposite mass is given by a zigzag (armchair) boundary 
whose properties can be related to those of the infinite 
graphene zigzag (armchair) ribbon with inhomogenous 
mass. For that matter, the top panel of Figs. [7] and [8] 
shows the current density map obtained for the infinite 
ribbon with At = — A^. Whereas a priori this current 
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FIG. 8: (Color online) (a) Current density map corresponding 
to the armchair infinite ribbon (along the interface, see text), 
(b), (c) and (d) as in Fig. (tI for the case of zigzag ribbons 
with opposite mass A_l ~ ^Ar — 0.3 t, under a magnetic 
field with $ = 0.0081 {B = 640 T). 



density should not be identical to the boundary of Fig. 
fl|a), due to the finite width of the of the Hall bar, our 
results indicate that it permits to anticipate the existence 
of backscattering induced by interface states. 

Since we are interested in the interface states as a 
source of backscattering, it is not necessary to consider a 
long central region with disorder, as we did in the previ- 
ous section. In addition, the role of disorder here could be 
to reduce the efficiency of the interface states to produce 
backscattering, improving the conductance thereby. 

For both armchair and zigzag geometries the conduc- 
tance is zero for £' < A — 0.3 t and is quantized in 
G = ^ above the gap, over an energy interval that co- 
incides with the absence of interface states, shown in the 
panel (a) of Figs. [TJandls] The comparison of the current 
density for an infinite ribbon with the two terminal con- 
ductance for both heterostructures, reveals a relation be- 
tween the existence of interface states in the junction ((a) 
panels) and the backscattering in transport ((b) panels). 
In particular, the two spectral regions with null current 
density at the interface give quantized plateaus of con- 
ductance. This is particularly apparent in the case the 
plateau with G — -f^ ior E right above the band gap. 

The connection between interface states and backscat- 
tering if further confirmed by studying the conductance 
as a function of the width of the Hall bar, at four differ- 
ent energies. In Figs. ItFc) andlslc) at the energies black 
and blue, the backscattering is completely canceled for 



sufficiently wide ribbons. At those energies, there are no 
interface states. In contrast, for energies red and green, 
the conductance oscillates as a function of the ribbon 
width. This can be interpreted as follows. As the rib- 
bon width is increased, the discrete spectrum of interface 
states shifts. When a interface state is in resonance with 
the electrode states, the backscattering is possible and 
conductance is reduced. In contrast, peaks in the trans- 
mission corresponds to poor matching between incoming 
state and interface state. It is apparent that the ampli- 
tude of the oscillations does not decrease significantly as 
the width of the ribbon increases, even in the presence 
of disorder. This suggest that the localization length of 
the interface states along the direction perpendicular to 
transport is longer than the ribbon width. The results 
are qualitatively similar in the case of a heterostructure 
made of zigzag ribbon with opposite masses. In this case 
the domain wall separating the two regions with opposite 
mass runs along the armchair direction. 

Disorder seems to have two effects on the interface- 
state induced backscattering in these heterojunction. 
One one side, is probably increasing the mixing of edge 
states to interface states, which should enhance the 
backscattering. On the other side, for sufficiently wide 
ribbons disorder could result in the localization of the 
interface states that are responsible of backscattering, 
which should decrease the backscattering. The compar- 
ison of the two curves with and without disorder in the 
middle panels of Figs. [7] and [8] indicate that, the domi- 
nance of one effect over the other depends on energy. In 
general, the interface induced backscattering effect is not 
qualitatively affected by disorder. 



VI. DISCUSSION AND CONCLUSIONS 

We have studied the electronic properties of graphene 
quantum Hall bars with a position dependent mass A. 
We have considered the case of stepwise constant A. We 
have found that at the boundary of two regions with op- 
posite A both in-gap kink states and interface states ap- 
pear. Interface states arise from the mixing of counter- 
propagating edge states that coexist in energy. In con- 
trast, kink states arise in the domain wall between two 
gapped regions and do not coexist in energy with bulk 
states. In the case of zigzag ribbons, there is one kink 
state at each valley, whose propagation direction changes 
from valley to valley. 

We have studied transport in two different configura- 
tions that would permit to probe either the in-gap kink 
states or the interface states. The study of transport 
parallel to the interface between two graphene regions 
with opposite A in the quantum Hall regime, permits 
one to study the properties of the kink states (section 
III). We have found that the coexistence in real space 
of two counter-propagating kink states, corresponding to 
the two valleys, leaves them unprotected from backscat- 
tering created by disorder. In the case of zigzag ribbon, 
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backscattering requires changing valley, which in turn re- 
quires short-range scattering, provided by the Anderson 
disorder. 

The study of transport across the junction of two semi- 
infinite ribbons with opposite mass, discussed in sec- 
tion V, permits studying the effect of interface states as 
sources of backscattering. Our calculations show how an 
incoming electron to the junction from a chiral edge state 
could scatter to a kink state propagating from one edge 
the opposite, enabling backscattering at the specific en- 
ergies at which interface states exist. 

Our calculations represent a toy model for situations 
in which graphene quantum Hall bars have a position de- 
pendent mass. This could be the case of a mass driven by 
electronic order, for which different ground states could 
coexist in the sample, or a mass modulated by the inter- 
action with a substrate with a very large commensuration 
period.1^1133 
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other. A central quantity in the method is the Green 
function G{EI — H) = I, where / is the identity matrix. 
The projection of the Green function operator over the 
central region can be written, after some algebra, as: 



Gc 



EI - H, 



c 



(A4) 



where the self energies S,, of the rj = L,R lead are given 
by 



^r, == V^g^V^ 



(A5) 



and gjj = ( -E/jj — H,j 1 are the projections of the Green 

function operators over the rj — L, R spaces. 

The conductance can be calculated in the linear re- 
sponse regime, within the Landauer formalism as a func- 
tion of the energy E. In terms of the Green function of 
the system, 33 it reads 



2p2 2e2 

G = '-^T (E) ^ '-^Tr 



^lGc^rG^j 



(A6) 



where T (E), is the transmission function across the con- 
ductor, and r^ — i[S,, — EJj] is the coupling between the 
conductor and the rj — L, R lead. 



Determination of the electrode Green function 



Appendix A: Calculation of Transmission 

1. Partition Method and Green functions 

In this appendix we give the technical details of the 
calculation of the Transmission function T{E), which 
yields the two terminal conductance through the Lan- 
dauer formula.*'^ The calculation method applies for a 
one dimensional system that can be split in three re- 
gions, a central "device" of finite size described with the 
Hamiltonian matrix H^ coupled to two semi-infinite elec- 
trodes, left and right leads, described with Hl and Hr. 
In matrix form the Hamiltonian reads: 



The Eq.(A4) to Eq.(A6l are all expressed in terms of 



the electrode Green function g^. In particular, when rep- 
resented in a local basis, it is the so called surface term 
of the Qri matrices which is needed. When the electrode 
Hamiltonian is written in the form: 



„t 



H,, 







V. 



V 



\ 



/ 



(A7) 



it can be shown that the electrode Green function satisfies 
the self-consistent equation: 



where 



and 



H 



Hs 



He V 

Ft Hs 



Hl 
Hr 



(Al) 



(A2) 



gr,= {E-hri- Vr^gr^vl) 



(A8) 



and corresponds to the central and one of the most time 
consuming steps in the calculation. We label the elec- 
trode Green function obtained in step i of the iteration 
procedure, as g* . We have found that the stability of the 
self-consistent procedure is improved by using the follow- 
ing algorithm to compute the step i, for i > 1 



V 



Vr V, 



R 



(A3) 



with Vl and Vr the coupling to the left and right leads, 
which we assume to be otherwise decoupled from each 



9r) 



a5r'+/35r' + (i 



■P)9. 



i-3 



(A9) 



with the the initial guess 5° 
/3 are mixing parameters. 



g and a. 
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3. Calculation of the transmission 

The other source of computational overhead in the cal- 
culation of Eq. ( A6 ) is the inversion of the central region 



He matrix, renormalized with the self-energies, to ob- 
tain Gc- However, this can be greatly simplified taking 
advantage of two facts. First, in the computation of the 
transmission, only a few matrix elements of the device 
Green function are actually needed, in particular, those 



involved in the F^Gc products, which are a minor frac- 
tion, given the surface nature of the F matrices. Second, 
the device Hamiltonian can be written as a tridiagonal 
block matrix. This permits to use specific techniques 
for tridiagonal matrices that make the procedure much 
faster.ESl Taking advantage of this approach, it is possi- 
ble to compute the transmission of 300 nm long bars in 
a desktop computer. 
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